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Abstract 

We examine the evolution of a universe comprising two interacting 
fluids, which interact via a term proportional to the product of their 
densities. In the case of two matter fluids it is shown that the ratio of 
the densities tends to a constant after an initial cooling-off period. We 
then obtain a complete solution for the cosmological constant (w = — 1) 
scenario. Finally, we investigate the general case in which the dark energy 
equation of state is p = tup, where w is a constant, and show that periodic 
solutions can occur if w < —1. We further demonstrate that the ratio of 
the dark matter to dark energy densities is confined to a bounded interval, 
and that this ratio can be O(l) at infinitely many times in the history of 
the universe, thus solving the coincidence problem. 

PACS: 98.80.-k, 95.36.+X 

I. INTRODUCTION 

Observations based on Type la supernovae and the cosmic microwave back- 
ground suggest that the universe consists mainly of a non-gravitating type of 
matter called 'dark energy', as well as a substantial amount of gravitating non- 
baryonic 'dark matter' [1-5]. Surprisingly, although the density of dark matter 
is expected to decrease at a faster rate than the density of dark energy through- 
out the history of the universe, their magnitudes are comparable today. This is 
known as the 'coincidence problem', and various attempts at its solution include 
the use of tracker fields [6] and oscillating dark energy models [7] . 

We will discuss a third possibility that has gained some attention recently, 
which is that dark energy and dark matter interact via an additional coupling 
term in the fluid equations [8-69]. The interaction is usually assumed to take 
the form AHp m + BHpA, where H is the Hubble parameter, p m and p\ are the 
dark matter and dark energy densities respectively, and A, B are dimensionless 
constants. However, in this paper we will instead consider an interaction of the 



* S.Z.W.Lip@damtp.cam.ac.uk 



1 



form 7p m /?A, where 7 is a (non-dimcnsionless) constant [8]. Such an interaction 
is natural and physically viable, since we would expect the interaction rate to 
vanish if one of the densities is zero, and to increase with each of the densi- 
ties. This form of interaction has also been used to model systems ranging from 
standard two-body chemical reactions to predator-prey systems in biology. Sta- 
tistical fits to observed data suggest that this form of coupling helps to alleviate 
the coincidence problem [9] , and it has been shown that among holographic dark 
energy models with an interaction term ^p m p K (where a, f3 ^ are integers), 
the one with a = /3 = 1 gives the best fit to observations [10]. 

We propose to investigate this model in more detail. Our analysis will con- 
sider the interaction of a dust fluid and a second fluid with an equation of state 
p = wp, where w is a constant. If the second fluid is also dust (w = 0), we show 
that the ratio of the densities tends to a constant after an initial cooling-off pe- 
riod. Thus, matter fluids coupled in this way can be considered, at late times, 
to evolve as a single non-self-interacting matter fluid. We then go on to obtain 
a complete solution of the system for conventional (p = — p) dark energy, but 
show that such a model cannot address the coincidence problem. Finally, we 
exhibit the various scenarios that can arise for other forms of dark energy, and 
demonstrate that if w < — 1 and there is an energy transfer from dark energy to 
dark matter, periodic solutions can arise in which the dark densities are com- 
parable for a substantial fraction of the evolution of the universe. This would 
address the coincidence problem. 

II. SETTING THE STAGE 

We work in a spatially flat Friedmann-Robertson- Walker universe, and as- 
sume that it contains the following perfect fluids: dark matter, with density 
p m , pressure p m and equation of state p m — 0; 'dark energy', with density p\, 
pressure pa and an equation of state of the form pa = wpAC 2 , where w is a con- 
stant (in the case w = — 1, we have a standard cosmological constant); baryonic 
matter, with density pb, pressure pb and equation of state pb — 0; and radiation, 
with density p r , pressure p r and equation of state p r — p r c 2 /3. 

The conservation equations governing the evolution of these fluids are: 



Pm = -3Hp m + 7p TO pA (1) 

PA = -3(1 + w)HpA - -fprnPA (2) 

Pb = -'SHpb (3) 

p r = -4Hp r (4) 

3H 2 = 8irG(p m + pa + Pb + Pr) (5) 



where H — a /a is the Hubble parameter, and a(t) is the scale-factor. Here, an 
overdot indicates a derivative with respect to (cosmic) time t. Note that the 
constant 7 is not dimensionlcss; it has dimensions of volume per unit mass per 
unit time. If 7 > 0, energy is transferred from dark energy to dark matter; the 
opposite occurs if 7 < 0. In the case 7 = 0, the two fluids do not interact. In [11] , 
it is pointed out that 7 < would worsen the coincidence problem, and in [12], 
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it is argued that we need 7 > in order for the second law of thermodynamics 
and Le Chatelier's principle to hold. These results are borne out later in Section 
5 of this paper. Note that we do not include baryonic matter in the interaction, 
due to the constraints imposed by local gravity measurements [70-71]. 

From the above equations, it is clear that pb and p r decrease as a~ 3 and a~ A 
respectively. In the next three sections of this paper, we will restrict ourselves 
to a late-time analysis, where pb and p r are small and can be neglected. This 
leads to the reduced set of equations: 



p rn = -'iHp m + JpmPA (6) 

PA = -3(1 + w)HpA - 'JpmPA (7) 

3H 2 = 8irG( Pm + p A ) (8) 

For simplicity, we will also adopt units in which 8nG — 1 and c = 1 in these 



sections, which contain a qualitative, mathematical analysis of the system. We 
will return to the full set of equations in Section VI when attempting to constrain 
the parameters of the system using observations. 

Differentiating ((5J) with respect to t and substituting for p m and pa gives 
the auxiliary equation 

2H = -p m - (l+w)p A (9) 

which will be useful later. Note that this equation is independent of the coupling 
parameter 7; this is a consequence of energy conservation. 

It will be of interest to consider the acceleration of the universe, which is 
given by the relation a = a(H + H 2 ). From ([8]) and ([9]), we have 

S = -§0»m + (l + 3tt;)pA). (10) 


Hence, since we assume that a(t), p m and pa are always non- negative, it is 
necessary (but not sufficient) that w < — -| in order for the universe to accelerate 
today. 

III. A SIMPLE CASE: TWO DUST FLUIDS 

We start by considering a model in which our fluids are both dust. The 
equations of the system are: 

pi = -3Hp 1 +jp 1 p 2 (U) 

P2 = -IP1P1 - ZHp 2 (12) 

3H 2 = Pl +p 2 (13) 

Define u = pi+ p 2 - Then ii = —^/3u 3 ^ 2 , and integrating gives ^= = y3(t + C), 
where C is an arbitrary constant. Therefore, u = 3 ( t + c y 1 and by shifting the 
time coordinate (t — > t — C) we can set C = 0. Hence H = ^= = ^ (as would 
be expected for a matter-dominated universe). 
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Substituting this result into (|TT|) and (fT3")) yields the equation 

2 / 4 



Pi = -7P1 +7Pi ( ^72 - Pi ) (14) 



which integrates to p x = t 2 (4Ae 4 7 /3t +3) , P2 = t 2( 4Be -4 7 /3t +3) ■ The condition 
Pi + P2 = 372- leads to the requirement that 16AB — 9, which leads to the 
general solution 

4 4 
Pi - 7-77757 , on'P2 



t2( 4 ^ e 4 7 /3t + 3) ' r* t 2 ^ e -4 7 /3t + 3) ■ 
The ratio of these densities is therefore 

Pi 3 _4-y 

- — = — e 3t 
p 2 

and, as t — > 00, this tends exponentially quickly to the constant 3/4A, with the 
individual densities evolving as pi ~ i~ 2 . 

This suggests that even if there were several types of 'dust' in the universe 
mutually interacting in this manner, we could approximate the evolution of their 
densities by treating them as non-interacting after a short initial period. (In 
other words, we might write pi ss Cit~ 2 for t > t* , where the relative magnitudes 
of the Ci are determined from a relatively short initial evolution up to time t* .) 

IV. 'CONVENTIONAL' DARK ENERGY (to = -1) 

We now consider a universe containing a dark matter fluid and a dark energy 
fluid, in which the dark energy behaves like a cosmological constant (that is, it 
satisfies the equation of state p\ = —p\). The dynamical equations simplify to: 

pm = -3Hp m + JPmPA (15) 
PA = -JPmPA (16) 

H 2 = \{ Pm + PK) (17) 

2H = -p m (18) 

In the non-interacting case 7 = 0, it is easy to see that the dark energy den- 
sity pa stays constant, and p m < (in fact, p m oc a~ 3 oc t" 2 ), so the dark 
matter becomes more diffuse. The ratio — decreases as t -2 , and so the coinci- 
dence problem remains. Since p m — > 0, we can see from (|10[) that the universe 
accelerates at an increasing rate as time progresses. 

Now suppose 7^0. First, observe that H = implies p m = p\ = 0, so we 
restrict attention to positive H. This implies an eternally expanding universe. 
From (HU) and (18]) it follows that p A = yle 27 ^ , where A is a positive constant. 
Given 7, we can find A based on the values of p\ and H observed today. Using 
(I17p . we can then parametrize the whole system in terms of H: 

Pm = 3H 2 -Ae 2 ~< H (19) 
PA = Ae 2 '< H (20) 



4 



From P0|) we can see that p\ never vanishes, so dark energy is a perpetual com- 
ponent of the universe. Now, using (1181) . we can obtain a first-order separable 
differential equation for H , and its solution is 




2dH 
Ae 2 ~< H - 3iJ 2 ' 



Unfortunately, the integral on the right hand side cannot be integrated using 
analytic methods. 

We now examine the acceleration of the universe in the two cases depending 
on the sign of 7. We can see that H + H 2 — —\H 2 + ^e 2 ~< H , so the universe 
accelerates iff Ae 2 ^ H > H 2 . In the case 7 < 0, we have p m < and pA > 0. 
Since the acceleration takes the same sign as 2p\ — p m , the acceleration occurs 
as a faster rate as time proceeds, just as in the non-interacting case. 

In the case 7 > 0, p\ < always, but the behaviour of p rn is more com- 
plicated. A representative example with 7 = 1 is given in Figure 1; for other 
values of 7, the qualitative behaviour is similar (in the sense that p m generally 
increases to a maximum and then decreases again). The straight line in the 
diagram corresponds to 2p\ = p m ; in the region of phase space above this line 
the universe is accelerating, and in the region below this line it is decelerating. 
It appears that in this model the universe can only go through at most a single 
decelerating phase. 
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Figure 1: Phase-plane diagram showing the evolution of the dark matter 
density p m and the dark energy density p\, for an interaction term p m pA 
(7 = 1). In this case, the dark energy is assumed to behave like a 
cosmological constant, and obeys the equation of state p = —p. The 
straight line represents universes with zero acceleration. The region of 
phase space above the line corresponds to an accelerating universe, and 
the region below corresponds to a decelerating universe. Note that this 
figure (and the next four) are intended only to display the qualitative 
behaviour of the model, so the units on the axes are arbitrary. 



From the expressions above we can see that r = = j[*}-,h — 1, and this 
quantity lies in the interval [—1, A j $ — 1]. (This can be seen by considering 

o 2 

the graph of y = — 1.) Thus, in order for the model to be feasible at all, 

we require the constants to satisfy the constraint Aj 2 < 3e~ 2 . 

In order for this model to address the coincidence problem, we require that 
r ~ O(l), and thus H 2 and e 2lH should be of comparable magnitude. We argue, 
however, that this cannot happen. 

In the case 7 < 0, it is easy to see that the ratio r always decreases with 
time, because p m ^ and p\ ^ at all times. Consider, then, the case 7 > 0. 
Suppose p m > c ^ at all times. Then H will decrease at a rate of at least | 
per unit time (by (fTS)) ). and will reach eventually - but this corresponds to 
pm = 0, which is a contradiction. Therefore p m must approach asymptotically 
at late times, so that the quantity 3H 2 — Ae 2lH becomes arbitrarily small, and 
comparable to p\ = Ae 2jH (which must therefore also be arbitrarily small). 
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Therefore 3H 2 , and therefore H, becomes arbitrarily small. This is impossible 
because 3H 2 - Ae 2 ~< H -> -A as H -> 0. 

Hence, this model cannot alleviate the coincidence problem. 

V. GENERAL DARK ENERGY 

Now suppose that the dark energy has equation of state p = wp, where w 
is a constant, and define K = w + 1. The previous case corresponds to K = 0, 
and the case in Section 3 (two dust fluids) can be regarded as the case K = 1. 
The general interacting case is then described by the following equations: 

p m = ~3Hp m + JPmPA (21) 
Pa = -JpmPA - 3KHp\ (22) 

3H 2 = p m + PA (23) 

We assume that 7^0. Adding ([21]) and (J22J) yields 

Pm + PA = -\/3(p m + Pa)(Pto + Kpa) (24) 
and subtracting appropriate multiples of (|2"Tj) and (|2"2"j) yields 

PmKpA ~ jKp m PA = PAPm + IplnPA- (25) 



Combining (|24[) and (|25|) . and integrating, gives the following relation between 
the dark energy and dark matter densities: 

pA = p K e2j V^TpI/Vs_ (26) 

Due to this rather awkward relation between the two densities, it is difficult 
to make further progress using purely analytical methods. However, we can 
still examine the cosmic evolution by treating (f2"Tj) - (f2"3|) as a two-dimensional 
dynamical system: 

Pm = -V^VPm + PAPm + IPmPA 
PA = -IPmPA - VS\/Pm + PaK p\. 

We consider the fixed points of this dynamical system. The values of p m 
and p\ at these fixed points must satisfy p m {lPA — \/3(p m + pa)) — and 
pAilPm + Ky/3(p m + pa)) = 0. It is clear that the origin (p m ,pA) = (0,0) is 
always a fixed point. For fixed points other than the origin, we need the other 
two factors to both vanish, and this leads to the consideration of the following 
four cases based on the signs of 7 and K: 

A. Case 1: K > 0,7 > 

There is a single fixed point (0,0). We can see that pa ^ always, so the 
dark energy density never increases. If pa is sufficiently great to begin with, the 
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matter density increases and then decreases. The condition p rn = corresponds 

2 2 

to the parabolic nullcline p m = 3 — p\. 

An illustration of the evolution is shown in Figure 2. In this, as well as the 
other cases, the actual values of 7 and K do not seem to affect the qualitative 
behaviour of the orbits, although their signs do. 

This scenario is unlikely to represent our present universe. Although for 
most trajectories there is an initial dark energy dominated phase which gives 
way to a dark-matter dominated phase, the model does not exhibit the late dark 
energy dominated phase which we are currently experiencing. 
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Figure 2: As in Figure 1, with the same interaction term p m ph (7 = 1), 
but now assuming a dark energy component whose equation of state is 
P=-p/2. 



B. Case 2: K > 0,7 < 

There is a single fixed point (0,0). We can see that p m ^ always, so the 
dark matter density never increases. If p m is sufficiently great to begin with, 

the dark energy density increases and then decreases. The condition pj\ = 

2 2 

corresponds to the parabolic nullcline p^ — 7 „m ~ Pm ■ 
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In this scenario, an initial matter-dominated phase gives way to a dark en- 
ergy dominated phase, and then the density of dark energy falls steeply. How- 
ever, this model does not address the coincidence problem, because the ratio 
r = — keeps decreasing (if ^ 1). 




C. Case 3: K < 0,7 < 



In this situation, dark matter transfers energy to 'phantom' dark energy. 
As usual, there is a single fixed point (0,0). Since we always have p m ^ 
and pa 0, matter is continually being converted into dark energy, and p m 
tends monotonically to (as p\ tends monotonically to 00). As the ratio of the 
densities ^ is always decreasing, this aggravates the coincidence problem and 
substantiates the results from [11] and [12] alluded to earlier in Section 2. 
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Figure 4: As in Figure 1, but with an interaction term —p m PA (7 
— 1) and a phantom dark energy component with equation of state p 
-3p/2. 



D. Case 4: K < 0,7 > 

In this case, 'phantom' dark energy transfers energy to dark matter. The 
system now has two fixed points at 

(nn , , - v (3K{K-1) 3(1 -K) 

(0,0) and (p m , p A ) = I 
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The eigenvalues of the Jacobian at the second fixed point are ±2yK{l — K)\ 
they are both purely imaginary, and so {pm,ph) is a center. Phase-plane plots 
indicate the existence of stable periodic orbits around this center, with no con- 
vergence to nor divergence from (p m ,p A ). The two nullclines are the parabola 

22 22 
Pm = ^-3^ ~ Pm ( on which p m = 0) and the parabola p A = \j^2 — Pm (on which 

A linearization about the centre (p m ,PA) gives the equations 

3K 3KC2K-1) 

x = + Y 1 y 



3(K 
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where we have defined x = p m — p m , y = p\ — pa and made the assumption that 
x, y are both small. Fitting the equation of a conic section Ax 2 + By 2 + xy = 
const and differentiating gives a solution A = — ^, B = K — i. In order for 
the trajectories to be ellipses, we require that 4AB — 1 > 0, but this is true (as 

expected) since AAB -1 = 2 (y^K + -jL= ^ . 

Thus, for all K < 0, the orbits of the system are bounded. For a given 
trajectory, the ratio ^ will oscillate between two extremes given by the gradi- 
ents of the tangents from the origin to this trajectory. Note that the interval 
defined by these two extremes contains the value = which is 0(1) if K 
is not too close to 0. In such a model the dark energy and dark matter will be 
comparable at infinitely many times, and this would provide a solution to the 
coincidence problem. An example of the evolution is shown in Figure 5 for the 
parameter values 7 = 1 and K = 




2 4 6 8 

Pm 



Figure 5: As in Figure 1, with an interaction term p m PA (7 = 1) but 
a phantom dark energy component with equation of state p = —3p/2. 
Note that the system now exhibits periodic orbits, and that the ratio of 
the dark densities for any particular orbit is constrained to lie within a 
bounded interval. 

VI. CONSTRAINING THE MODEL PARAMETERS 

In the previous two sections, we have seen that, in the class of interacting 
models with an interaction of the form 7p m pA, the only model that could pos- 
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sibly help to alleviate the late-time coincidence problem is the one in Section 
V.D. Thus, we will now consider this model in more detail. 

The analysis in the previous sections holds only at late times, when the 
baryon and radiation densities are small and can be neglected. In order to 
constrain the parameters 7 and K of the model using observations of the past 
universe, we need to examine its past evolution. It is therefore necessary to 
include baryons and radiation in the analysis, since their densities are significant 
at early times. 

Before we proceed, we should emphasize that the following discussion is not 
an attempt to obtain the best-ht parameters, but merely to show that the model 
under consideration is a viable fit to observations for non-zero 7 and K. We 
have seen that the strength of this model is that it addresses the coincidence 
problem (which the standard concordance model cannot do), and we will find 
that it appears to fit observations at least as well as the concordance model. 
In order to obtain best-fit parameters and make a more detailed comparison to 
observations, more careful work is needed, including a full stability analysis of 
perturbations. This is, however, outside the scope of the present paper, and is 
left as a topic for future investigation. 

The complete system of equations to be solved is: 



where we have transformed the independent coordinate from cosmic time t to 
rcdshift z, and a prime indicates a derivative with respect to z. The quantities 
Pbo and p r0 are the present densities of baryons and photons, respectively. 

In order to find suitable choices of parameters for 7 and K 7 we can integrate 
these equations numerically for various values of these parameters, and compare 
the results to observations. In particular, we shall attempt to examine the effect 
of three observational constraints which are claimed to be model-independent 
[72] : (i) the shift parameter R, related to the angular scale of the first acoustic 
peak in the CMB power spectrum, (ii) the distance parameter A, related to the 
measurement of the BAO peak from a sample of SDSS luminous red galaxies, 
and (iii) Type IA Supernovae data. 

The shift parameter R is defined as follows [72, 73]: 



where z r ccomb = 1091.3 [74] (WMAP7) is the recombination redshift. When 
testing our model, we will adopt the values n m o = 0.227 for the present frac- 
tional density of dark matter, fi^o = 0.0456 for the present fractional density 
of baryons, and Ho — 70. 4km s _1 Mpc -1 for the present value of the Hubble 
parameter. Note that this is only an approximation, because the derivation of 



(1 + z)p' r , 



(1 + z)Pa 
3H 2 



rn 




8?rG(p m + pa + p6o(l + zf + p r0 (l + zf) 
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these values from the WMAP7 data assumes a standard, non-interacting ACDM 
model. The value of R obtained from the WMAP7 data is 1.725 ± 0.018. 
The distance parameter A is defined as follows [75, 76]: 



where z b = 0.35. The value of A has been determined to be 0.469(n s /0.98)"°' 35 ± 
0.017 (la constraints) [76], and we will take the scalar spectral index n s to be 
0.963 [74] (WMAP7). We will also assume the WMAP7 parameters for O m0 
and flbo (discussed above). 

Using these parameters, we have calculated R and A for various values of 
7 > and K < 0. The observational constraints from R put a crude upper 
bound on 7 of about 1.4 x 10 7 m 3 kg -1 s _1 , and requires K > —0.48. The 
constraint from A suggests that K > —0.27 and requires 7 to be less than 
about 1.8 x 10 9 m 3 kg -1 s _1 . The closer K is to 0, the higher the values of 
7 permitted by either test: a more negative value of K tends to increase the 
values of both A and R. 

In the remainder of this paper, we will take the parameters 7 = 100 and K = 
—0.15. Note that, although these are not necessarily the best-fit parameters, 
they are in accord with the observed values of R and A. 

We have fitted this particular model to SNIa data [77]. At the outset, we 
might expect that there will not be a large deviation from the behaviour of the 
concordance model in this regime, since the densities of either dark matter or 
dark energy would be very low at small rcdshifts, and the interaction term would 
therefore be negligible. The recent evolution of the system would therefore be 
similar to that of a standard ACDM universe. The results confirm this: for a 
sample of 608 Union2 supernovae, the model with 7 — 100, K = —0.15 gives a \ 2 
value of 756.4, which is slightly better than the one obtained for the 7 = 0, K = 
model (780.7). 



We now consider the concrete instance of our model with 7 = 100 and 
K = —0.15. The evolution of the fractional densities with redshift is shown in 
Figure 6: 




VII. DISCUSSION 
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Figure 6: The past evolution of the fractional densities Qi with redshift, for a 
universe in which 7 = 100 and K = —0.15. The universe undergoes the usual 
sequence of radiation, dark matter and dark energy dominated eras. 

From Figure 6, we can see that this model provides a suitable sequence 
of cosmological eras corresponding to what we know about our universe. The 
universe is radiation-dominated at very early times. This is followed by a long 
matter-dominated era, with fl\ <C Q m , up to around redshift z « 0.5. At late 
times, we encounter a period of dark-energy-dominated acceleration. 

The past history of the universe is somewhat dependent on the parameter 7. 
In particular, increasing 7 to larger than about 10 induces a peak in the baryon 
fractional density at about redshift 2500 and leads to a brief baryon-dominated 
era around that time, whilst moving the transition from the radiation-dominated 
era to the matter-dominated era later in time (i.e., to smaller redshift). 

Also, note the peak in the fractional density of dark energy at around z 
2500. This corresponds to the beginning of a periodic cycle. There are, in 
fact, infinitely many such cycles, as the evolution continues into the past, and 
the oscillations get more rapid as we get closer and closer to the Big Bang 
(since, at early times, the Hubble parameter is large due to the high density of 
radiation). However, this is not easy to see even if we extrapolate Figure 7 to 
higher redshifts, since at high redshifts p m and p\ are swamped by the radiation 
density. We can, however, observe this effect on a logarithmic plot (Figure 7): 
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Figure 7: A log-log plot of the fractional densities fii against redshift, for a 
universe with 7 = 100 and K = —0.15. This plot shows that many cycles may 
occur before the current one. 

Having established that the model adequately describes the past universe, 
let us consider its implications for the late universe. At late times, the densities 
of the non-interacting baryons and photons will decrease and become negligible. 
The evolution of the universe thus converges to a stable limit cycle, in which 
the following four phases repeatedly occur: 

• Dark energy is converted quickly into dark matter, while H remains large. 

• When the density of dark energy is sufficiently low, the interaction term 
in the evolution equation for p m is overwhelmed by the —2>Hp m term, 
and the density of dark matter starts to decrease, while the dark energy 
density stays low. The value of H decreases, and the universe decelerates. 

• The dark matter and dark energy densities are small and comparable. The 
dark energy density slowly starts to increase while the dark matter density 
continues to decrease. This marks the transition from a dark-matter- 
dominated universe to a dark-energy-dominated universe, and the universe 
begins to accelerate. The universe that we are currently experiencing is 
in the later stages of this transition. 

• When the density of dark energy then becomes sufficiently large, it con- 
tinues to increase until the interaction term becomes non-negligible. 
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This process goes on forever, with eternally alternating dark-matter- and dark- 
energy-dominated eras. Also, the universe undergoes infinitely many periods of 
acceleration and deceleration, as can be seen from equation ([TU)) . 

In this limit cycle, the dark matter and dark energy densities are comparable 
for a significant proportion of the time. The following graph shows a plot 
demonstrating this: 




Figure 8: The future evolution of the fractional densities tti with cosmic time 
t, for 7 = 100 and K = — 0.15. In this plot, t is measured in seconds, and the 
time t = corresponds to the present. 

The duration of an entire cycle is dependent on the choice of parameters: 
choosing a more negative value of K leads to more and shorter cycles. From 
Figure 8, the densities are comparable for about a quarter of each cycle, and 
that dark energy is dominant at most other times. This suggests that, in this 
model, it is not unnatural for the universe to be in a state at which the dark 
matter and dark energies are comparable. 

This model, however, has a limitation: it does not explain the apparent 
coincidence of the baryon and dark matter densities being comparable at the 
present time. In particular, the baryon density continues to decrease steadily 
as the cycles progress, whereas the minimum value of the dark matter density 
stays at roughly the same level over all future cycles, because of its periodic 
interaction with dark energy. 

VIII. CONCLUSIONS 
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We have considered a model with an interaction term proportional to the 
product of the densities of the interacting fluids, and applied it to the interaction 
of a dust fluid and another fluid with equation of state p = wp, for various values 
of the constant w. In the case of two dust fluids, we have obtained a general 
solution, and have shown that their density evolution is similar to that in the 
non-interacting case, after an initial short cooling-off period. 

We then proceeded to show that such an interaction between dark matter 
and dark energy in the form of a cosmological constant would not solve the 
coincidence problem. In the case 7 < 0, energy is transferred from dark matter 
to dark energy, and the ratio of the energy densities r = p m / PA always de- 
creases with time. We would therefore expect the dark matter density today 
to be negligible compared to the dark energy density. However, this does not 
correspond to what we observe in the universe today. On the other hand, in 
the case 7 > where energy is transferred from dark energy to dark matter, 
we found that it was still not possible to have the densities be of comparable 
magnitude. This led to the conclusion that such an interaction would not help 
to solve the coincidence problem. The same is true if 7 < and w is allowed to 
vary arbitrarily, or if w > — 1 and 7 is allowed to vary arbitrarily. 

However, it is interesting to note that if the dark energy follows a 'phantom' 
equation of state (w < —1) and energy is transferred from dark energy to 
dark matter, it is possible to obtain periodic orbit solutions. These correspond 
to cyclic situations in which the ratio of the dark densities is comparable at 
infinitely many times. We have also shown that suitable parameters 7 and K 
can be chosen so that the model is consistent with observations. 

Such a scenario could alleviate the coincidence problem, because it can be 
shown that if w is not too close to —1, the periodic orbits would enclose a fixed 
point corresponding to a density ratio r that is 0(1), and thus the value of r 
on these trajectories would be 0(1) at infinitely many times in the evolution of 
the universe. 
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